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ABSTRACT 

We  study  various  stability  type  conditions  on  a matrix  A related  to 
existence  of  a positive  diagonal  matrix  D such  that  the  Lyapunov  matrix 
O = AD  + DA*"  is  positive  definite.  Such  problems  arise  in  mathematical 
economics,  in  the  study  of  time-invariant  continuous-time  systems  and  in  the 
study  of  predator-prey  systems.  Using  a theorem  of  the  alternative,  a 
characterization  is  given  for  all  A such  that  a corresponding  D exists.  In 
addition,  some  necessary  conditions  for  consistency  and  some  related  ideas 
are  discussed.  Finally,  a method  for  constructing  D is  given  for  matrices 
A satisfying  certain  conditions. 
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AH  the  vectors  in  this  paper  are  real  and  the  matrices  are  square 
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and  real. 
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The  well  known  theorem  of  Lyapunov  [ 1892]  states  that  all  the 


eigenvalues  of  a matrix  A have  negative  (positive)  real  parts  if  and  | 
only  if  there  exists  a symmetric  positive  definite  matrix  H such  that] 


0 = AH  + HA 

is  negative  (positive)  definite.  Such  a matrix  A is  said  to  be  negative 


A 


(positive)  stable  and  such  matrices  have  been  studied  extensively  in  the 
economics  and  mathematics  literatures  (e.g.  Johnson  [ 1974]  and  Taussky 
[ I96l  ]).  By  a stable  matrix  we  shall  mean  a positive  stable  matrix. 

We  are  concerned  here  with  conditions  under  which  the  matrix  H in 
the  expression  for  Q can  be  chosen  to  be  a positive  diagonal  matrix  D. 
When  this  is  the  case,  A is  called  diagonally  stable. 
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Matrices  A for  which  such  a matrix  D exists  are  important  in  the  study 
of  dynamical  systems  (e.g.  Araki  [ 1975]),  in  the  study  of  predator-prey 
systems  (e.g.  Krikorian  [ 1976 ])  and  in  the  study  of  dynamic  equilibra 
in  economics  (e.g.  Johnson  [l974]). 

We  first  collect  some  observations. 

Lemma  j.  The  following  statements  are  equivalent  for  a matrix  A. 

(l)  There  exists  a positive  diagonal  matrix  D such  that 

0 = AD  + DA^ 

is  positive  definite. 

da^  There  exists  a positive  diagonal  matrix  D such  that 

x^AD  X > 0 

for  ail  nonzero  vectors  x. 

(1^)  There  exists  a positive  diagonal  matrix  D such  that  the 
symmetric  part  of 

1 1 
D"  *AD* 

is  positive  definite. 

(l^)  A^  satisfies  Condition  (i). 

(1^)  A ^ exists  and  satisfies  Condition  (i). 

(Ig)  For  each  positive  diagonal  matrix  E,  AE  and  EA  satisfy 
Condition  (1 ). 

(1^)  Every  principal  submatrix  of  A satisfies  Condition  ( l ). 
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To 


Proof.  Clearly  Conditions  (l)  and  (1  J are  equivalent. 

d 

complete  the  proof  it  is  enough  to  show  that  (l)  implies  (l^)» 

(1^)  and  (1  p. 

Let  D be  a positive  diagonal  matrix  such  that  AD  + DA  = P 
is  positive  definite.  Then 

d”^ad^  + (d"^adV  = d“^(ad  + da'^)d"^  = d"^pd'2 

which  is  positive  definite,  so  that  (1  ) holds.  Clearly  (1  ) follows  since 

D C 

P^  is  positive  definite.  Also  A is  nonsingular,  being  stable,  and 

a"^d  + d(a“^)^  = a“^p(a"M^ 

which  is  positive  definite  so  that  (l^)  holds.  If  E is  a positive  diagonal 
matrix  then 

(AE)(e'^D)  + (e"^D)(AE)^  = P , 

(EA)(DE)  + (DE)(EA)^  = EPE 

so  that  (1^)  holds.  The  statement  (l)  (ip  follows  from  the  fact  that 
ptUsdpal  *4ibmatrlc»s  of  a positive  definite  matrix  are  also  positive  definite. 

We  now  characterize  these  matrices  utilizing  the  theory  of  the 
alternative  for  systems  of  linear  inequalities. 

Theorem  1.  A matrix  A is  diagonally  stable  if  and  only  if  for  every 
nonzero  positive  semi  definite  matrix  B,  BA  has  a positive  diagonal  element. 

Proof.  Let  be  the  space  of  real  diagonal  matrices  ordered  by 

the  cone  K of  nonnegative  diagonal  matrices.  Let  V 
1 ^ 


be  the  space  of  symmetric  matrices  ordered  by  the  cone  of  positive  semi- 
definite  matrices.  In  both  cases  we  use  the  inner  product 

(A,  B)  = tr(AB)  . 

If  the  dual  cone  K of  K C V,  is  defined  by 

K*  = {B  e V;(B,  A)  > 0 for  all  A e K} 
then  and  (cf.  Hall  [1967]).  Define  the  linear 

operator  ^2 

T^(D)  = AD  + Da'^  . 

The  theorem  of  the  alternative  (Berman  and  Ben-Israel  [1971],  Berman  [1973], 
p.  20)  says  that  the  system 

T^(D)  e int  K^,  D e int 

is  consistent  if  and  only  if  there  is  no  nonzero  solution  for  the  system 

-T^(B)  e Kj*,  B e . 

The  mapping  T^  is  defined  by 

(T^(B),D)  = (B,T^(D)) 

for  all  D e V . Thus  T^(B)  is  a diagonal  matrix  which  satisfies 
1 A 

tr(D^(B)D)  = tr(BT^(D))  = tr(B(AD  + DA^)  = tr((«A  + A^B)D) 
for  every  diagonal  matrix  D.  Thus 

T^(B)  = diag  (BA  + A^  B),^  = 2 diag  (BA).^  . 

The  conclusion  follows.  ■ 
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Similar  results  can  be  obtained  for  the  operators  defined  by  Stein 
[ 19  52),  Schneider  [ 1965  ] and  Hill  (1969).  For  the  Stein  operator 
we  have 

CAD)  = D - A*'  DA  and  C^(B)  = diagi  B. , - (ABA^  ), , } . 

A A 11  ' ji 

We  now  list  some  immediate  corollaries  to  Theorem  1.  The  first 
corollary  and  the  "only  if"  part  of  the  second  corollary  can  also  be 
found  in  Johnson  [1974]. 

Corollary  1.  If  A is  a triangular  matrix  with  positive  diagonal 
elements,  then  A satisfies  (1). 

Proof.  By  Lemma  1,  A can  be  assumed  to  be  upper  triangular.  Let 
k be  the  smallest  index  such  that  for  a given  positive  semidefinite 
matrix  B,  b^^j^  > 0.  Then  (BA)j^j^  >0.  ■ 

Corollary  2.  A matrix  A of  order  2 is  diagonally  stable  if  and  only 
if  all  its  principal  minors  are  positive. 

Proof.  The  "only  if"  part  holds  for  matrices  of  any  order  as  we 
shall  show  in  the  next  theorem.  It  is  enough  to  prove  the  "if"  statement 
for  matrices  A of  the  form 
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We  must  show  that  for  every  nonzero  positive  semidefinite  B,  either 
(BA)  > 0 or  (BA)  > 0.  If  b..  = 0,  then  (BA).,>0.  If  b >0 

11  Cct  11  ii  Ij 

and  b^2  ^ ^ i^Qy  assume  that 


B = 


b 

c 


2 

with  0 < b < c.  We  also  have 
If  now 


1 > a 


12^21’ 


whence  c > 


(BA)ii  = 1 +a2^b<0  and 
(BA)22  = c + a^  2^  < 0,  then 
cK-a^jbK-a^^b)  = , 


a contradiction.  ■ 


Some  fairly  obvious  necessary  conditions  for  a matrix  to  satisfy 

\ 

Condition  (l)  are  now  given. 


Theorem  2.  If  a matrix  A satisfies  Condition  (l)  of  Lemma  1 then 
it  also  satisfies  each  of  the  following  conditions. 

(2)  All  the  principal  submatrices  of  A are  stable. 

(2  ) The  real  part  of  each  eigenvalue  of  each  principal  submatrix 

a 

of  A is  positive. 

(3)  All  the  principal  minors  of  A are  positive. 

( 3 ) Every  real  eigenvalue  of  each  principal  submatrix  of  A is 

a 

positive. 

( A + E is  nonsingular  for  each  nonnegative  diagonal  matrix  E. 
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(3  ) For  each  x ^ 0 there  is  a positive  diagonal  matrix  E such  that 

x^AEx  > 0 . 

(3  ,)  A does  not  reverse  the  sign  of  any  vector.  That  is,  if  x * 0 
d 

then  for  some  subscript  i 

x (Ax),  > 0 . 
i 1 

(3  ) For  each  signature  matrix  S = diag(±  l),  there  exists  an 
e 

X > 0 such  that 

SASx  > 0 . 

Proof.  It  is  well  known  that  Conditions  (2)  and  (2  ) are  equivalent 

■ 3 

by  the  theorem  of  Lyapunov  [1892]  for  arbitrary  matrices  A.  Moreover 
Conditions  (3)  through  (3^)  are  also  equivalent  (cf.  Plemmons  [ 1976]) 
and  (2  ) implies  (3).  Also,  it  is  clear  that  (l)  implies  (2)  by  Condition  (1  ).  ■ 

For  the  next  result  we  need  the  following  notation.  Let 
<n)  = {l,...,n}.  For  any  nxn  matrix  C and  any  aC<n>,  let 
C[<t>]  = 1,  and  for  a t let  C[a]  be  the  principal  submatrix  of  C 
with  indices  from  a and  let  C(a)  be  the  principal  submatrix  of  C 
with  indices  from  (n)  - a.  If  D = C ^ then  from  Sylvester's  identity, 

(cf.  Gantmacher  and  Krein  [I960],  p.  15)  we  have 

(*)  det  D[a]  = (det  C)"^  det  C(o)  . 

Theorem  3.  If  a matrix  A satisfies  any  of  Conditions  (1),  (2)  or  (3) 
in  Theorem  2,  then  A is  nonsingular  and  each  of  the  matrices  A^  and 
A ^ satisfies  the  same  condition. 
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Proof.  The  only  nonimmediate  assertion  is  that  if  A satisfies  (Z) 


then  so  does  A . Suppose  A satisfies  (Z).  Let  z be  any  complex 
number  such  that  Re  z ^ 0 and  put 

A = A + zl,  B = a"^  . 
z z z 

Clearly  every  principal  submatrix  of  A is  stable.  Then  for  any  a C (n) 
we  have  from  {*)  that 

det  B [a]  = (det  A ) ^ det  A (a)  . 
z z z 

But  the  stability  of  A (a)  yields  that  the  right  hand  side  is  nonzero. 

z 

Since  z is  any  number  in  the  closed  right  half  plane  it  follows  that 
B[a]  has  no  eigenvalues  in  the  closed  left  half  plane,  i.e.  B satisfies  (2).  ■ 
For  a given  A we  have  the  following  result  on  the  inductive  construc- 
tion of  a D satisfying  (1). 


Theorem  4 . Let  A be  partitioned  (n  - 1, 1)  so  that 

/A, 


A = 


T 


\ a 


Suppose  there  is  a positive  diagonal  matrix  D,  of  order  n - 1 such  that 
AjDj  + DjA|  = is  positive  definite.  If  z = (z^)  is  a vector  of  order 
n - 1,  define  the  inner  product  (2^  1^2)  by 


Let  X = a^,  y = D^a^.  Then  a positive  scalar  d exists  such  that  for 
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D = 


Di  0 


, AD  + DA  is  positive  definite,  if  and  only  if 
a > ((xly>  + (x|x>  (y  (y>)  . 


AD  4 DA  = 


Proof.  Regarding  d as  an  unknown  we  note  that 

a^d  + D,aJ 

is  positive  definite  if  and  only 


AjDi  4 D^Aj  a^d  + D^a^ 
2ad 


a‘Di  +da‘ 


if  its  determinant  A is  positive.  But  A = (det  A^^)[(2da  - a^D  + da^) 


^ ^1^3^  ^ ^ ^ when  and  only  when 
0 < 2da  - - da^^A^'^D^a 


3 ■ '*^2*11  ®2'* 


= 2da  - (yly>  - 2<ylx)d  - <xlx>d^  . 

2 

Put  p(d)  = (x|x)d  4-  (2(x|y)  - 2a)d  + (yly).  There  is  a d for  which 
p(d)  <0  if  and  only  if  for  the  discrlmant  of  p we  have 

0 < 4{((x  |y>  - a)^  - <x  |x>  <y  ly)  } . 

2 2 

Put  q(z)  = z - 2<x|y)  z 4- (<x|y)  - (x|x)  <y|y>).  We  have 
q'(z)  = 2z  - 2<x  |y>,  so  at  z - (x|y),  q(z)  has  the  negative  minimum 
q(<x  |y>)  = -<x  lx)  (yly).  The  larger  root  of  q(z)  = 0 is  z = <x|y>4- 
(x  lx)  (y  ly)  . So  for  a>a^  = (x|y)  + (x|x)(y|y)  the  discriminant  of  p 
is  positive.  In  order  that  a d > 0 exists  we  must  have  that  the  larger 
root  of  p(d)  = 0 is  positive.  However  the  larger  root  of  p(d)  = 0 is 
given  by 

r(a)  = a - (x|y)  4-  ((x|y)  - a)^  - <x|x)  (y|y)  . 
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We  see  that  r(a  ) > 0 and  r(a)  is  increasing  for  a > a . Hence  for 
o o 

any  a > a^  the  discriminant  of  p is  positive  and  the  larger  root  is 
positive.  Thus  we  can  find  the  required  d.  ■ 

Remarks. 

1.  If  A satisfies  any  of  Conditions  (1),  {^)  or  (3)  then  so  does 
each  principal  submatrix  of  A. 

If  A satisfies  Condition  (I)  of  Lemma  1,  then  there  need  not 
always  exist  a positive  diagonal  matrix  D such  that 

AD  + DA^  = I . 

In  fact,  the  proof  of  Theorem  3 in  Taussky  [1961]  shows  that  such  a D 
exists  if  and  only  if 

A = E + SE 

where  E is  a positive  diagonal  matrix  and  S is  skew  symmetric. 

3.  It  is  easy  to  construct  matrices  which  satisfy  (3)  but  not  (2) 
of  Theorem  2 and  thus  not  (1).  For  example,  the  matrix 


satisfies  (3),  but  its  eigenvalues  are  9 and  ±'V3i  so  that  it  is  not 
stable  and  thus  does  not  satisfy  (2).  Note  also  that  A is  normal. 

4.  To  show  that  (2)  does  not  imply  (1)  we  proceed  as  follows:  Let 
A be  the  matrix  given  in  Remark  3.  Let 
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with  t > 0.  Then  for  sufficiently  small  <,  all  the  principal  submatrices  of 
AD  are  stable.  Thus  AD  satisfies  Condition  (2),  but  by  Remark  3 and 
Lemma  1,  AD  does  not  satisfy  Condition  (l). 

5.  The  matrices  A = (a_)  arising  in  the  work  of  Krikorian  [ 1976], 

mentioned  earlier,  satisfy  >0  for  all  i * j.  The  example  in 

Remark  4 shows  that  even  for  these  matrices  Condition  (2)  does  not  imply 
Condition  (1). 

6.  For  certain  classes  of  matrices.  Conditions  (1),  (2)  and  (3) 
are  all  equivalent.  These  classes  include: 

(a)  The  triangular  matrices. 

(b)  Ihe  positive  definite  matrices. 

(c)  The  2X2  matrices. 

(d)  The  matrices  A=  (a^^)  with  < 0 for  i # j. 

The  matrices  satisfying  (d)  and  (1),  (2)  or  (3)  are  of  course  the  non- 
singular M-matrices. 

7.  A significantly  large  class  of  matrices  which  satisfy  (1)  can 


be  obtained  as  follows:  For  A 

= (a„), 

let 

???(A)  be  the  matrix  given  by 

^ 

r 

if 

i = j 

Ij  I 

L-K.I 

if 

1 * ) 
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Then  it  is  known  (e.g.  Plemmons  [1976])  that  if  A has  a positive 
diagonal  then  each  of  the  following  conditions  is  eguivalent: 

(a)  ^(A)  is  a nonsingular  M-matrix. 

(b)  There  exists  x > 0 such  that  for  each  signature  matrix 
S = diag(±  I ) , 

SASx  > 0 . 

(c)  There  exists  a positive  diagonal  matrix  D such  that  AD  is 
strictly  row  diagonally  dominant. 

(d)  There  exists  a positive  diagonal  matrix  D such  that 

D~^AD 

is  strictly  row  diagonally  dominant. 

It  follows  from  the  work  of  Johnson  [ 1974]  that  any  matrix  A 

with  a positive  diagonal,  satisfying  any  of  these  conditions,  also 

satisfies  (1).  That  the  converse  is  not  true  is  easily  established. 

Let 


Then  by  Corollary  2,  A satisfies  (1).  But  V}{A)  is  a singular  matrix. 
Open  Question 

A matrix  A was  defined  by  Arrow  and  McManus  [ 1958]  to  be 
D-stable  if  A satisfies: 
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(4)  For  each  D = diag(d^,  . . AD  is  stable  if  and  only  if 

d.  > 0,  i = 1, . . n. 

The  more  common  practice  now  is  to  define  A to  be  D-stable  if  it 
satisfies 


( S)  AD  is  stable  for  each  positive  diagonal  matrix  D. 

The  class  of  matrices  satisfying  ( S)  properly  contains  the  class  of 
matrices  satisfying  (4)  which,  in  turn,  properly  contains  the  class 
satisfying  (1).  That  these  classes  are  not  the  same  is  illustrated  by  the 
matrix 


A = 


(o 


-1 


which  satisfies  (4),  but  not  (1).  Examples  also  exist  that  satisfy  (5) 
but  not  (4)  (e.g.  Johnson  [1974]). 

In  view  of  these  conditions  we  pose  the  following  open  question: 
Are  the  topological  interiors  of  the  sets  of  matrices  satisfying 


(1),  (4)  and  (5)  identical? 
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